Adsorption and dissociation of H 2 monomer on ceria(lll): Density functional 

theory calculations 
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The adsorption properties of isolated H2O molecule on stoichiometric and reduced ceria(llll) 
surfaces are theoretically investigated by first-principles calculations and molecular dynamics simu- 
lations. We find that the most stable adsorption configurations form two hydrogen bonds between 
the adsorbate and substrate. The water molecule is very inert on the stoichiometric surface unless up 
to a high temperature of 600 K. For the reduced surface, we find that the oxygen vacancy enhances 
the interaction. Moreover, simulations at low temperature 100 K confirm that it is facilitated for 
water to dissociate into H and OH species. 

PACS numbers: 68.43.Bc, 68.43.Fg, 68.43.Jk, 68.47.Gh 



I. INTRODUCTION 

Ceria (Ce02) is a crucial material for catalytic applica- 
tions due to its remarkable properties 1 ^. One of its im- 
portant application is that ceria can effectively promote 
the water-gas shift reaction, CO + H2O — > CO2 + H2, 
by which converting harmful carbon monoxide to the 
less harmful carbon dioxide, and generating hydrogen. 
A full account of the involved mechanism has not been 
achieved, while the interactions between water and ceria 
surface may plan a key role during the reaction. 

Therefore, a systematic study of the fundamental in- 
teractions involved in the water-ceria system is much de- 
sirable. The oxygen-terminated ceria(lll) surface, which 
is most often found in experiments 3 -, has been identified 
to be energetically the most stable 4 ^— . The formation of 
oxygen vacancies in ceria, hence reducing the cerium from 
Ce 4+ to Ce 3+ , is closely related to its interesting proper- 
ties and important implications, such as oxygen storage 
capacity of ceria 1 -. A heavily debated question is on the 
favorable position of the oxygen vacancy. Theoretically, 
different research methods including density-functional 
theory (DFT) calculations, have given different favorable 
vacancy sites ( on-surface or on-subsurface)^. For the 
interaction between adsorbate and ceria surface, obvi- 
ously, the on-surface oxygen vacancy is much concerned. 

Experimentally, it has been found that water on re- 
duced ceria powders can effectively oxidize Ce 3+ into 
Ce 4+ , with the H2 production 1 ^—. More recently, evi- 
dences have been provided by Henderson et alJ£ for that 
the water may split at low temperature below 170 K, 
and exposure of water at 650 K results in additional sur- 
face reduction of ceria instead of oxidation. In addition, 
they offered three potential adsorption geometries for wa- 
ter on the stoichiometric cerial(lll) surface, namely, the 
so-called C2V geometry without hydrogen bond between 
the water and the substrate (no-H-band), single hydro- 



gen bond (one-H-bond), and two hydrogen bonds (two- 
H-bond) configurations, needing to be examined by the- 
oretical calculations. 

Several theoretical works have investigated the H2O- 
ceria interaction by DFT calculations. Kumar et al^- 
reported that the water adsorption on the stoichiometric 
ceria(lll) surface favors one-H-bond configuration, con- 
sistent with the later literature result o 15 i 16 . While this 
conflicts with the results obtained by Fronzi et alrH, who 
suggested a two-H-bond structure. For water on reduced 
(with on-surface oxygen vacancy) ceria(lll) surface, the 
calculated stable adsorption structures varied from no- 
H-bon d 16 ' 17 , to one-H-bond^ configuration, while Chen 
et al^ found no binding between water and reduced sur- 
face. Moreover, Kumar et alJ^ predicted that oxidation 
of the reduced surface by water with the production of 
hydrogen gas is weakly exothermic, while Yang et alJ£ 
and Fronzi et al/H argued that hydroxyl surface forms 
upon water dissociating. 

Such discrepancies may partially arise from the dif- 
ferent unit cell adopted by different authors. Actually, 
the water coverage adopted in previous works ranges 
0.25-1.0 ML. Too high water coverage would be prob- 
lematic because it may mask the true H20-ceria interac- 
tion, by the dipole-dipole interactions between adjacent 
water molecules. Therefore, a systematical study of wa- 
ter on stoichiometric and reduced ceria(lll) surfaces at 
lower coverage, which minimizes the lateral interaction 
between the adsorbates, is indispensable to shed a light 
on such discrepancies for a proper understanding of water 
adsorption on ceria(lll) surface. 

The DFT+ U calculations based on first-principles 
method have proved to be an effective approach for 
studying the strong correlation materials consisting d 
and / electrons, compared with the plain DFT method 18 . 
For bulk ceria, there has been a vast literature deal- 
ing with this issue, despite the controversial determi- 
nation of the Hubbard parameter U 6 i 14 ^ 19 ~ —. For the 



2 



surface properties of ceria, however, by declaring its neg- 
ligible influence, little work takes into consideration the 
U parameter— ~—. Nevertheless, it has been revealed that 
the plain DFT method failed to correctly describe the be- 
havior of the water molecule on reduced ceria surfac o 16 ' 17 . 
Due to the localization of / electrons, especially for re- 
duced ceria surfaces with oxygen vacancies, to describe 
more accurately the excess electrons, strong correlation 
effect must be considere d 14 ' 19 , hence the DFT+ U ap- 
proach is required. In addition to the static DFT calcu- 
lation, ab initio molecular dynamics (AIMD) simulations 
can effectively demonstrate the dynamical properties of 
a system when considering temperature, which can be 
compared with the existing experimental measurements. 
Moreover, the use of AIMD can be a beneficial comple- 
ment to the static calculations, especially for systems 
showing multiple low-lying minima in energy^. 

In the present work, we systematically study the prop- 
erties of isolated water adsorption on stoichiometric and 
reduced ceria(lll) surfaces at low coverage by means of 
first-principles methods. We perform the microscopic 
static DFT calculations aiming at elucidating the elec- 
tronic structure and static stability for the water adsorp- 
tion structures. Besides, AIMD simulations at various 
temperatures are also conducted to investigate the dy- 
namical properties of these systems. 



II. COMPUTATIONAL METHOD 

The calculations are performed within DFT using the 
projector-augmented wave (PAW) method, as imple- 
mented in the plane- Wave based Vienna ab-initio simu- 
lation package (VASP)^. The cutoff energy for the plane 
wave expansion is set to 400 eV. For the bulk ceria, the 
exchange and correlation effects are treated within both 
the local density approximation (LDA) and the PBE 24 
generalized gradient approximation (GGA). For the ceria 
surface, however, only LDA is adopted. Here the cerium 
5 s, 5p, 5d, 4/, and 6 s, and the oxygen 2 s and 2p electrons 
are treated as valence electrons. We use the DFT+ U for- 
malism formulated by Dudarev et al^ to account for the 
strong on-site Coulomb repulsion amongst the localized 
Ce 4/ electrons. In this scheme the on-site two-electron 
integrals, expressed in terms of two parameters, i.e., the 
Hubbard parameter U, which reflects the strength of the 
on-site Coulomb interaction, and the parameter J, which 
adjusts the strength of the exchange interaction, are com- 
bined into a single parameter U e ff = U — J, indicating 
that only the difference between U and J is significant. 
For the bulk ceria, the Coulomb U is treated as a vari- 
able, while the exchange energy J is set to be 0.7 eV. 
Then the parameter U e ff is fixed at a reasonable specific 
value for the surface calculations. 

The oxygen-terminated ceria(lll) surface is modeled 
by a slab composing of six atomic layers and a vacuum 
region of 15 A. A (3 x 3) surface unit cell, in which each 
monolayer contains nine atoms, is adopted in the study 



of the H2O adsorption. The isolated H2O molecule is 
placed on one side of the slab, equivalent to 1/9 cover- 
age. Our test calculations have showed that such low a 
coverage is sufficient to avoid the interaction between ad- 
jacent H2O molecules, hence makes the true H2 0-ceria 
interaction emerge. During our calculations, the bottom 
two atomic layers of the substrate are fixed, and other 
atoms as well as the H2O molecule are free to relax until 
the forces on the ions are less than 0.03 eV/A. Integration 
over the Brillouin zone is done using the Monkhorst-Pack 
scheme^. For bulk ceria, we use a 9 x 9 x 9 grid points, 
and for the surface we use a 3 x 3 x 1 grid points. 

The calculations of the energy barriers for the water 
diffusion and dissociation processes are performed us- 
ing the climbing image nudged elastic band (CI-NEB) 
method^ 7 ., which is an improved version of the tradi- 
tional NEB method^ for finding minimum energy re- 
action paths between two known minimum energy sites, 
resulting in a more accurate estimation of the activation 
energy than the regular NEB method does, by driving up 
one of the intermediate states near the transition point 
to reach the highest saddle point along the reaction path. 

For AIMD calculations, the Born-Oppenheimer ap- 
proximation is used; ions are considered as classical ob- 
jects moving on the potential surface created by the elec- 
trons obeying the quantum mechanical equation. The 
canonical ensemble using the Nose thermostat^ is em- 
ployed. Electronic energies are converged up to 10 -6 eV. 
With a time step of 1 fs, initial equilibration steps are 
performed over 1 ps, while production runs are 6 ps long. 



III. RESULTS AND DISCUSSION 
A. Determination of bulk properties of ceria 

The stoichiometric ceria crystallizes in a CaF2-like 
ionic structure with space group Fm3m. The experimen- 
tally determined lattice parameter is about 5.41 A 30 ' 31 , 
and the bulk modulus varies over a rather broad range of 
204 to 236 GPa22~— . In this section we demonstrate the 
dependence of the lattice parameter clq and bulk modu- 
lus B on U c ff in the range 0^8 eV, obtained by fitting 
the third-order Birch-Murnaghan equation of stated, as 
presented in Fig. [1] 

We can see from Fig. Q] that the reasonable lattice 
parameter values ao are 5.40 A (LDA+C/, U e ff=5 eV) 
and 5.48 A (GGA+ U, U cff =3 eV), and the bulk mod- 
ulus values B are 210.3 GPa (LDA+ U) and 178.4 GPa 
(GGA+J7). Noticeably, our calculations are consistent 
with previous studies of Andersson et al^ and Loschen 
et al^, but we can not reproduce the results obtained by 
Jiang et al&. It is clear that the GGA+ U calculations 
overestimate the lattice parameter and underestimated 
the bulk modulus, and thus give a slightly less accurate 
description of the atomic structure of ceria. On the other 
hand, LDA+ U proposes a reasonable description. More- 
over, within the latter choice, the calculated 02p-Ce4/ 
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FIG. 1: Dependence of the lattice parameter (a) and bulk 
modulus (b) of CeC>2 on U e ff. The dotted lines stand for the 
experimental measurements. 
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FIG. 2: (Color online) PDOS for the H 2 molecule and the 
first-layer O and the second-layer Ce of the stoichiometric 
surface for (a) free H2O and the clean ceria(lll) surface, (b) 
the Al adsorption site. The inset in (b) shows the 3D electron 
density difference, with the isosurface value set at ±0.05 e/A 3 . 
The area filled with yellow color represents MOs of H2O. The 
Fermi level is set to zero. 



and 02p-Ce5ci band gaps are 2.21 and 5.29 eV, con- 
sistent with the experimental measurements of about 3 
and 6 eV, respectively^. Moreover, it has been revealed 
that the LDA treatment is actually better than the GGA 
ones when simulating chemical reactions on reduced ceria 
surfaces 19 . In addition, as mentioned above, most of the 
previous DFT calculations dealt with the H^O-ceria in- 
teractions by using GGA functional^—, with somewhat 
conflicting results existing, hence the LDA treatment can 
be performed as a beneficial complement. Therefore, in 
the following surface calculations, we will present only 
the LDA+ U results. 



B. H2O adsorption on stoichiometric ceria(lll) 

The structural and energetic parameters of the free 
water molecule are calculated within a box with the same 
size of the adsorbed systems. The optimized geometry for 
free H2O gives a bond length of 0.97 A and a bond angle 
of 105.2°, consistent with the experimental values of 0.96 
A and 104.4°^. 

The adsorption energy of the system is calculated as 
follows: 

E a d = £'H 2 0/ceria(lll) ~ -E-H 2 — £-ccria(lll) > (1) 

where £h 2 o, ^ceria(ui), an d ^HsO/ccriaCm) are the to- 
tal energies of the H2O molecule, the clean ceria(lll) 
surface, and the adsorption system respectively. Accord- 
ing to this definition, a negative value of indicates 
that the adsorption is exothermic (stable) with respect 
to a free H2O molecule and a positive value indicates 
endothermic (unstable) reaction. 



Adsorption geometry optimizations are performed for 
a variety of initial orientations of the H2O molecule 
placed over different surface sites, i.e., the flat-lying (wa- 
ter molecule lying parallel to the surface), O-down (O 
atom pointing to the surface) and H-down (two H atoms 
pointing to the surface) structures. We identify different 
H2O adsorption structures, whereas only the low-energy 
one (Al) is discussed here, with the adsorption energy 
of —0.96 eV and the corresponding frequencies of 3269.7, 
3206.8, and 1490.3 cm -1 . The equilibrium geometry of 
Al state is illustrated in the inset of Fig. [2] in which 
the adsorbed H2O monomer lies nearly parallel to the 
surface, with the O atom of water located almost on 
top of the second-layer Ce atom, while H atoms sym- 
metrically oriented pointing to the first-layer O atoms 
(do-H=l-85 A), namely, with two-H-bond configuration 
formed, in agreement with the results of Fronzi et alM>. 
Noteworthily, this adsorption configuration is quite dif- 
ferent from previous studies, either the arbitrary assump- 
tion of the no-H-bondi^ or the one-H-boncU^— configu- 
ration. The difference may arise from that those previ- 
ous studies underestimated the hydrogen bond strength, 
which plays a key role in the interactions between water 
and oxide surfaces. Moreover, the dipole-dipole interac- 
tions between adsorbed water molecules at high coverage 
probably blur the exact monomer adsorption configura- 
tion studied in this letter. 

In order to further understand the precise nature of 
the chemisorbed molecular state, we plot in Fig. [2] 
the electronic projected density of states (PDOS) of the 
molecular adsorption structure for the stable adsorp- 
tion configuration of Al in comparison with those of 
free H2O molecule and clean ceria(lll) surface. The 
three-dimensional (3D) electron density difference Ap(r), 
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FIG. 3: (Color online) Variation in H2O adsorption energy on 
the stoichiometric ceria(lll) surface as a function of the lat- 
eral displacements of O atom from its original adsorption site 
Al to the nearest-neighbor one. Insets show the correspond- 
ing geometry sketches for minima (Al, A2) and transition 
(Tl, T2) states. The chart summarizes the targeted O-H dis- 
tances (in A) pointed out in the insets. Large, moderate and 
small spheres stand for Ce, O and H atoms, respectively. 



which is obtained by subtracting the electron densi- 
ties of noninteracting component systems, p coria ( m )(r) + 
/? H2 °(r), from the density p(r) of the H20/ceria(lll) sys- 
tem, while retaining the atomic positions of the compo- 
nent systems at the same location as in H20/cerai(lll), 
is also shown in the inset of Fig. [2] (b). Positive (blue) 
Ap(r) indicates accumulation of electron density upon 
binding, while a negative (yellow) one corresponds to 
electron density depletion. Molecular orbital (MO) 2a\ 
of water (not shown here) is far below the Fermi level 
and thus remains intact in water-metal interaction. Here 
we consider only three MOs I62, 3ai, and lb\. 

As illustrated in Fig. [3J upon adsorption, the orbital 
I62 is observed to shift down in energy by 0.47 eV. More 
significantly, the orbitals 3ai and l&i undergo noticeable 
changes, especially for the latter, becoming more broad- 
ened and hence apparently delocalized, with the domain 
of energy corresponding to that of the ceria surface. Ob- 
viously, the d- and /-states of the Ce atom do not change 
much because of localization. On the other hand, the 
shifting up in energy for the peak of 0-2p states of the 
ceria surface play a key role in the water-ceria interac- 
tion, which indicates that the adsorbed MO I61 may act 
as an electron donor state. The features of the orbital hy- 
bridization are further substantiated by the 3D electron 
density difference plotted in the inset of Fig. H^b). We 
can see that there exists a large charge accumulation be- 
neath the O atom of the adsorbate. This accumulation 
should come from the charge redistribution because of 
derealization of the MOs of water, rather than from the 
substrate, which is consistent with above PDOS analy- 
sis demonstrating the localization of Ce electronic states. 
Moreover, obvious charge transfer from H atoms of the 
adsorbate to the surface O atoms is observed, suggest- 
ing strong hydrogen-bond interactions between the ad- 
sorbate and substrate. 

Given a H2O molecule at a stable adsorption site, it is 



interesting to see how it diffuses on the substrate. There- 
fore, we calculate the diffusion path and energetic barrier 
of water on ceria(lll) surface between neighboring sta- 
ble Al sites, with the adsorption energy as a function of 
the lateral displacement of O atom of the water molecule 
shown in Fig. [3] As indicated, the molecular Al state 
is adsorbed by —0.96 eV. Along the diffusion path, the 
system evolves through an energy barrier of 0.35 eV to 
reach a semi-stable state A2, which is also nearly flat on 
the surface. A2 is less stable with an adsorption energy 
of —0.63 eV, with the O atom of water located on top of 
the third-layer O atom, while as in the Al state, two hy- 
drogen bonds formed between H atoms of the water and 
the surface-layer O atoms (<io-H = l-80 A). At the tran- 
sition state Tl, one H atom of the water points to the 
top O atom, with the distance e?o-H = l-56 A. Along this 
path, the diffusion from A2 to the neighboring Al is hin- 
dered by a barrier of 0.18 eV, across the transition state 
T2, which also forms single hydrogen bond. It can be 
seen from the transition states Tl and T2 that the single 
hydrogen bond formed between the water monomer and 
the oxide surface is unfavored in energy, which is different 
from the previous high-coverage studies^ 4 - but consistent 
with our finding that the two-hydrogen-bond interaction 
is favored for the monomer adsorption. 

According to the overall energetic profile obtained 
above, we calculate the diffusion coefficient at differ- 
ent values of temperature T of this system within 
the Arrhenius-type relation, D — vod 2 e~ 6 / kBT , where 
vo=3269.7 cm -1 is the vibrational frequency of Al state, 
ri=4.90 A is the displacement of the water during diffu- 
sion, and e=0.51 eV is the activation energy. For com- 
parison the diffusion coefficient at T=300 K and 600 K 
are calculated to be 7.18 x 10" 14 and 1.30 x 10~ 9 m 2 /s, 
respectively. It is clear that at 300 K, the coefficient is 
too small to facilitate the diffusion process. While at 600 
K, diffusion should be observed easily. 

This is confirmed by our AIMD simulations, starting 
from the adsorbed molecular site Al at 300 K and 600 
K, which are shown in Fig. 2J At room temperature (300 
K), it can be seen that the molecular state Al is recovered 
with only small oscillations. Apparently the H atom os- 
cillates more intensely than the O atom, with the longest 
<io-H=l-54 A in the water molecule. At higher tempera- 
ture of 600 K, the water molecule undergoes more consid- 
erable oscillations (with the longest do-H=l-63 A), and 
finally negotiates a diffusion to the neighboring adsorp- 
tion site Al, as identified in the above CI- NEB calcula- 
tions. 

Moreover, during AIMD simulations we find that at 
temperature as high as 600 K, neither desorbed nor dis- 
sociated water configuration is observed. This indicates 
that on one hand the water-oxide interaction is strong 
enough to stand for high temperature, on the other hand, 
the stoichiometric ceria surface is inactive towards H2O 
splitting. Actually, by static calculations, we find that 
the co-adsorption of H and OH species on the regular ce- 
ria(lll) surface is endothermic thus unstable. In a word, 
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FIG. 4: (Color online) Projection on the stoichiometric ce- 
ria(lll) surface of the position of the oxygen atom (blue) 
corresponding to the H2O from AIMD simulation at the tem- 
perature of 300 K (a) and 600 K (b), respectively. The black 
arrowhead in (b) denotes the diffusion orientation. The large 
red and black filled circles mark the equilibrium positions of 
the first-layer O atoms and the second-layer Ce atoms, respec- 
tively, while the smaller red ones represent the third-layer O 
atoms. 




FIG. 5: (Color online) PDOS for (a) free H 2 and the clean 
reduced ceria(lll) surface, (b) the A2' adsorption site. The 
isosurface value for the 3D electron density difference is set at 
±0.05 e/A 3 . The red open circle in the inset of (b) denotes 
the oxygen vacancy position. The Fermi level is set to zero. 



the stoichiometric ceria(lll) surface does not dissociate 
water, in good agreement with previous calculations^— . 
However, the existence of defects, especially the oxygen 
vacancies on the surface of ceria, may act as active sites 
for the dissociation of water. Hence, in the following sec- 
tion we will focus our attention on this topic. 

C. H2O adsorption on reduced ceria(lll) 

It has been proposed that for bulk ceria, the oxygen- 
vacancy formation process is facilitated by a simultane- 
ous condensation of two electrons into localized /-level 
traps on two cerium atoms 37 . The situation is similar for 
ceria surface^. We can see from Fig. that in the re- 
duced surface the Ce-4/ states become partially occupied, 
meanwhile other states of the surface remain almost in- 
tact, compared with that in stoichiometric surface. This 
observation may contribute much to the adsorption prop- 
erties of water on the reduced ceria surface. 

Similarly to the stoichiometric surface, we explore dif- 
ferent H2O adsorption structures and find that the most 
stable molecular adsorption state is A2 adjacent to the 
defect (labeled A2') with the adsorption energy of —1.09 
eV, while Al adjacent to the defect is the next most 
stable state (labeled Al') with the adsorption energy of 
—0.93 eV. As illustrated in the insets of Fig. [fU for these 
two stable molecular states, the distances between the 
H atoms of water and first- layer O atoms are 1.77 and 
1.82 A, respectively, both smaller than that on the sto- 
ichiometric surface. In A2' configuration the O atom of 
the adsorbed water is located on top of the third-layer O 
atom, with a lower adsorption energy compared with that 
on the stoichiometric surface. Apparently, as depicted by 
previous studies^—, the oxygen vacancy enhances the 



water-ceria interaction. 

Figure EJb) presents the PDOS of the molecular state 
A2'. One can see that the MOs I62 and 3ai keep local- 
ized, only with a down-shift in energy by 0.16 and 1.36 
eV, respectively, while l&i becomes totally delocalized, 
which gives the main contribution to the adsorption be- 
havior of water monomer. In addition, the other contri- 
bution to the adsorption behavior comes from the 0-2p 
states of the surface, and the water adsorption introduces 
new peak for 2p states of the surface O atom, aligning 
in energy with 3a 1. From the 3D electron density differ- 
ence, we can observe more obvious charge transfer from 
H atoms of the water to the surface O atoms, compared 
with that on stoichiometric surface, probably providing 
a reason why the water molecule adsorbs to the reduced 
surface more strongly. Resembling that on the stoichio- 
metric surface, the partially filled Ce 4/ state almost 
does not participate the interaction with water, which 
exactly reveals the importance of the Hubbard parame- 
ter U for correctly depicting its localization. Moreover, 
we notice that this parameter prominently dominates the 
bulk properties of ceria such as lattice constant, hence 
indirectly holds the interactions between water and sub- 
strate, particularly for the reduced ceria(lll) surface. 

To clarify how an oxygen vacancy affects the dynamical 
properties of the adsorbed water monomer, let us exam- 
ine the diffusion and dissociation of the H2O molecule on 
this reduced surface. We calculate the diffusion path of 
the water between the stable states A2' and Al' (see Fig. 
[BJ. It is clear that the diffusion from A2' to Al' needs 
to overcome an energy barrier of 0.44 eV, while the in- 
verted diffusion from Al' to A2' is more facilitated with a 
lower barrier of 0.27 eV. Noticeably, the diffusion energy 
barriers presented here are larger than that on the sto- 
ichiometric ceria surface, and this should be attributed 
to the rotation of H2O molecule during the diffusion pro- 



6 




Reaction coordinate (A) 



FIG. 6: (Color online) Variation in energy for the diffusion 
and dissociation of H2O molecule on the reduced ceria(lll) 
surface. 
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FIG. 7: (Color online) Projection on the reduced ceria(lll) 
surface of the position of the oxygen atom corresponding to 
the OH species (purple) and the hydrogen atom (orange) from 
AIMD simulation at the temperature of 100 K (a) and 600 
K (b), respectively. The red open circle denotes the oxygen 
vacancy position. 



cess, which can be seen from the transition state T3 that 
has an O-down structure. 

To determine the final configurations of the dissocia- 
tion, we performed AIMD simulations for water on the 
reduced ceria(lll) surface, starting from the adsorbed 
molecular site A2', at different temperatures (see Fig. [7]). 
What is interesting, even at temperature as low as 100 K, 
is that the water monomer first oscillates around its equi- 
librium position and soon starts to dissociate into H and 
OH species almost in situ. It is found that the O atom of 
the OH species occupies the position of oxygen vacancy, 
with the O-H bond almost perpendicular to the surface, 
oscillating together with the surface. Whereas, another 
oscillating H atom of the dissociated water molecule is 
bonded to the adjacent top-layer O atom. Namely, hy- 
droxyl surface forms upon the dissociation of water on 
the reduced ceria surface. This is the same at higher 
temperature of 600 K, only with a more intensive oscil- 
lation. 

With this final dissociation configuration, now we re- 
turn to the static CI-NEB calculations. As presented in 



Fig. [6l the adsorption energy of the dissociated water 
state Dl was calculated to be —2.52 eV, with the bond 
length between H and O of 0.98 A. To reach the state Dl, 
the system only needs to overcome a small barrier, 0.02 
eV, corresponding to the low temperature 100K, as de- 
picted in the AIMD simulations. At the transition state 
T4 leading to the dissociation, the interatomic distance 
between the O atom of the OH species and another H 
atom is quite small, 1.12 A. These results consist well 
with previous experimental observation^ and theoreti- 
cal calculation^ 6 -, while differ from the calculation with 
plain DFTil that resulted in a high activation energy 
2.35 eV for water dissociation. 

Interestingly, by static calculations we have also iden- 
tified a semi-stable adsorption configuration with the ad- 
sorption energy —0.89 eV, in which the O atom of the 
water occupies the position of the oxygen vacancy, with 
a no-H-bond configuration formed. This structure has 
once been considered as the most stable one on the re- 
duced ceria surfac e 16 ' 17 , and also a specious initial state 
for evolving hydrogen gas^ at high coverage (0.5-1.0 
ML). However, we find that the process of H2 molecule 
production is energetically unfavorable when compared 
with the formation of hydroxyl surface (about 0.52 eV 
higher in energy), therefore it cannot be observed in the 
AIMD simulations. Furthermore, it has been pointed out 
that the energy needed for H2 production is as high as 
2.05 eV— , consistent with our results and the experi- 
mental observation that the reduced ceria(lll) surface is 
resistant to being oxidized by water—. 



IV. CONCLUSIONS 

In conclusion, we have systematically studied the ad- 
sorption and dissociation behaviors of H2O monomer on 
ceria(lll) surface by first-principles calculations comple- 
mented with AIMD simulations. The most stable adsorp- 
tion configurations on both the stoichiometric and re- 
duced surfaces were identified, with the electronic struc- 
tures analyzed in detail. Blocked by a diffusion barrier 
of 0.51 eV, water on the stoichiometric surface is very 
immobile up to at room temperature, while at higher 
temperature of 600 K, the water molecule can freely dif- 
fuse on the surface. For the reduced surface, We have 
found that the reaction of water and ceria is enhanced by 
the on-surface oxygen vacancy. It is facilitated for water 
to dissociate into H and OH species, only hindered by 
a small barrier of 0.02 eV, which was confirmed by the 
molecular dynamics simulations at temperature as low 
as 100 K. Besides, the reduced ceria(lll) surface is resis- 
tant to being oxidized by water, in good accordance with 
the existing experimental results. The present results 
obtained at low water coverage afford to settle existing 
discrepancies about the behaviors of water on ceria sur- 
face, and provide a guiding line for deeply understanding 
the water-ceria surface interactions. 
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